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Abstract 

We propose a mean field (MF) theory for a tiomogeneously driven granular gas of 
inelastic particles with Coulomb friction. The model contains three parameters, a 
normal restitution coefficient r„, a maximum tangential restitution coefficient r™, 
and a Coulomb friction coefficient The parameters can be tuned to explore a 
wide range of physical situations. In particular, the model contains the frequently 
used // — > oo limit as a special case. The MF theory is compared with the numerical 
simulations of a randomly driven monolayer of spheres for a wide range of parameter 
values. If the system is far away from the clustering instability (r„ ~ 1), we obtain a 
good agreement between mean field and simulations for /i = 0.5 and rj" = 0.4, but 
for much smaller values of r„ the agreement is less good. We discuss the reasons of 
this discrepancy and possible refinements of our computational scheme. 

Key words: Kinetic and transport theory of gases. Computational methods in 
fluid dynamics PACS: 47.50+d, 51.10.+y, 47.11. +j 



1 Introduction 

Granular gases [1] are usually described as collections of macroscopic parti- 
cles with rough surfaces and dissipative interactions. In order to study them, 
kinetic theories [2-4] and numerical simulations [5] were applied for special 
boundary conditions. The dynamics of the system is assumed to be dominated 
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by two-particle collisions, modeled by their asymptotic states: A collision is 
characterized by the velocities before and after the contact, and the contact 
is assumed to be instantaneous. In the simplest model, one describes inelastic 
collisions by normal restitution r„ only. However, surface roughness is im- 
portant [3,5], since it allows for an exchange of translational and rotational 
energy. Here we briefly sketch a study of a model where a Coulomb friction 
law with coefficient /i and a tangential restitution coefficient account for 
tangential inelasticity and friction [5,6]. We first introduce the model, then we 
describe a MF theory for the driven granular gas of rough spheres, and finally, 
we compare analytical results and numerical simulations. In the conclusions 
we discuss possible future refinements of the computational scheme. 



2 The Model 

We consider 3-dimensional spheres of mass m and diameter 2a interacting 
via a hard-core potential, confined on a 2-dimensional (2D) layer of linear size 
L, with periodic boundary conditions. Inelasticity and roughness are described 
by a normal restitution r„, a Coulomb friction law with friction /i, and a tan- 
gential restitution rt which depends on r^. A* the collision angle jc for 
sliding contacts and on a maximum tangential restitution for sticking con- 
tacts, when the tangential elasticity becomes important. When two particles 
1 and 2 collide, their velocities after collision depend on the velocities before 
collision through a collision matrix whose elements depend on r„, n, 7c, and 
rj". Thus we calculate the momentum change using a model that is consistent 
with experimental measurements [6]. From the momentum conservation laws 
for linear and angular direction, energy conservation, and Coulomb's law we 
get the change of linear momentum of particle 1 as a function of r„, fi, and 
[5]. The change of the normal component of the relative velocity depends on 
the normal restitution r„, which is a tunable parameter, while the change of 
the tangential component of the relative velocity depends on by the tangential 
restitution coefficient rt — min rf,r^ , where is the coefficient of maxi- 
mum tangential restitution, —1 <r^<l.The quantity rf is determined using 
Coulomb's law such that for solid spheres rf = —1 — (7/2)yu(l + rn) cot ■jc with 
the collision angle 7r/2 < 7^ < tt [5]. Here, we simplified the tangential contacts 
in the sense that exclusively Coulomb-type interactions, i.e. AP^*) is limited 
by ijlAP^"\ or sticking contacts with the maximum tangential restitution rf* 
are allowed [5]. 



3 The mean field Theory 

We start from the results of Huthmann and Zippelius [3] for a freely cooling 
gas of infinitely rough particles. Here we apply their MF theory to the case of a 



2 



driven gas of rough particles, which is the most common experimental situation 
[7]. The MF kinetic theory of Huthmann and Zippelius is formulated for a gas 
of rough particles with constant tangential restitution r^, corresponding to the 
limit fj, — oom our model. It is based on a pseudo-Liouville-operator formalism 
and on the assumption of a homogeneous state, with a Gaussian probability 
distribution of translational and rotational energies. The main outcome of 
this approach is a set of coupled evolution equations for the translational and 
rotational temperatures Ttr and T^ot [3]. Here, we write down and solve the 
MF equations for Ttj. and Tj-ot for a granular gas of rough particles in which 
the translational velocities are subjected to a random uncorrelated Gaussian 
driving of variance ^q- These equations read for a 2D layer of spheres 
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where r] — q{l + rt)/(2q + 2), q = 2/5 for spheres, n is the gas density and x is 
the pair correlation function at contact. 

We can use the Verlet-Levesque [8] approximation in 2D x= (1 ~ 70/16)/(l — 
0)^ where (j) is the volume fraction of the gas. For long times the system 
approaches a steady state. By imposing -^Tf^ — 0, ^T^ot = we get the equi- 
librium temperatures 

- m ( and ^ - - ^ (2) 

where Qd = 2%^/^ /T{D /2) = 2% ioi D = 2 and 7 = ^ + ^(l+rt) (6 - n)- 
(5/49)(l +r()^/(9 — 5rt) for spherical particles. By linearizing the set of Eqs. 
(1) around T^^"^ and T,^'^^ wc get the final approach to the steady state: 6Trot{t) — 
RSTtr{t) and Ttr{t) —Tl'^' = 5Ttr{t) ~ 5Tjr(0) exp[— 370;^]. The quantity uj — 

Q/jXna Enskog collision frequency for elastic particles at the 

temperature T^^, and tc — (370;)"^ oc 7"^/^ is a characteristic relaxation time 
[4]. The MF Eqs. (1) can be applied to the three parameter model in the limit 

^ = 00, see Fig. 1. 

However, experimental measurements are well reproduced by this model only 
for finite /i. Here we investigate the possibility to describe the effect of finite 
friction /i by replacing in the Eqs. (1) with its average (r^) over the proba- 
bility density P[sin(7t)] of the normalized impact parameter b/{2a) = sin(7t), 
with the condition < rj". The used simplifying approach (7c ~ 7i) ignores 
the fluctuations of due to its dependence on the rotational degree of freedom, 
and we expect that it is correct only in some trivial limiting cases. Neverthe- 
less, this simple approximation allows us to realize that for rj" = 0.4, fj, — 
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0.5, r„ ~ 0.9, corresponding to many experiments, simulations fit well with 
the modified MF theory, as we will se below. If the molecular chaos hypoth- 
esis is vahd P[sin(7()] = 1, and we get (r^) = — + r„) ln^\/l + c^ + c) , 
with c= 1(1 + r™)/(/^(l + r„)). 



4 The simulations 



Here we compare numerical simulations of a randomly driven monolayer of 
spheres, performed by using an Event Driven (ED) algorithm [5], with the 
MF predictions (for details see [4,5,9,10]). Every simulation is equilibrated 
without driving with r„ = 1 and rj" = —1. Then inelasticity and driving 
arc switched on. We used a fixed volume fraction = 0.34, = 11025 and 
different values of r„, r™, /i. In Fig. la-d, the translational and rotational tem- 
peratures for fixed r„ = 0.95, fi = 10^ and different values of arc rescaled 
with the MF equilibrium temperatures and plotted versus the rescaled time 
t/tc (X ^7^/^. The high value of fi allows to decouple normal and tangential 
momentum (r^ = r^). The agreement is good, even for low positive values 
of rj". For negative r™, agreement with MF is observed for the translational 
temperature, while the rotational temperature shows deviations from scaling 
in the transient phase, although the equilibrium value fits well with MF, for 
rj" ~ —1. This is due to a failure of the approximation STrot{t) ~ RSTtr{t) 
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Fig. 1. (a)-(b) Rescaled translational and rotational temperatures Ttr{t)/Tf^ , 
Trot{t)/T^ot versus the rescaled time t7^/^, for r„ = 0.95, /x = 10^ and positive 
fT'i (c)-(d) the same as (a)-(b) but for negative rj". 
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Fig. 3. (a)-(b) Rescaled translational and rotational temperatures Ttr{t)/Tf^, 
Trot{t) /T^ot versus the rescaled time ^7^^^, for r„ = 0.95, = 0.4 and different 
jjL] (c)-(d) the same as (a)-(b) but for ji = 0.5, = 0.4 and different r„. 

for weak coupling between rotational and translational degrees of freedom 
and finite time. In Fig. 2a-b we plot the same quantities as in Fig. 1 but for 
r„ = 0.95, = 0.5. We obtain for r™ near to unity significant deviations from 
mean field theory, while for negative r™ the agreement with mean field theory 
is very good. The explanation for this results is that for rj" ~ — 1 and high 
enough /x, one has (r^) ~ r[", and this correspond to have — > oo. In Fig. 3a-b 
we show simulation results for r„ = 0.95 and variable for = 0.4. These 
simulations confirm the previous interpretation. Finally, in Fig. 3c-d we show 
simulation results for ii = 0.5, = 0.4 and different r„. In this case the data 
collapse is very good although deviations from MF are observed. This is a 
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coincidence, since (rt) ~ 0.363 ~ rj". To improve the MF theory for the finite 
H case, the effect of a random must be fully taken into account. We actually 
are studying the possibility to include a collision angle dependent in the 
ensemble average of the energy variation due to collisions [3] . The MF theory 
we obtain seems to be very promising and gives a good qualitative agreement 
with simulations. A paper is in preparation with the results of this study [10] 
and also a three-dimensional analysis is in progress [11]. 
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